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Abstract 

We study low-temperature non Gaussian thermal fluctuations of a system of 
classical particles around a (hypothetical) crystalline ground state. These thermal 
fluctuations are described by the behaviour of a system of long range interacting 
charged dipoles at high-temperature and high-density. For the case of uniformly 
bounded fluctuations, the low-temperature linked cluster expansion describing the 
contribution to the free energy is derived and analysed. Finally some nonpertubative 
results on the existence and independence of boundary conditions of the Gibbs states 
for the associated dipole systems are obtained. 
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Chapter 1 
Introduction 



The question why and how crystals form at low temperature belong to the major open 
problems in theoretical physics [1,2,18]. 

Considering a system of N classical particles described by some Hamiltonian function 
H, we would like to know that the ground state configurations, which correspond to 
the local minima of H, are very regular and that under a small perturbation obtained 
by coupling this A^-particle system to a heat bath at temperature T^, the typical particle 
configurations for small values of do not fluctuate much away from the crystallic ground 
state conflguration. It is clear that certain surface effects could destroy this naive picture 
and therefore one has to pass to the thermodynamic limit which corresponds to letting 
N ^ oo while keeping the density of particles fixed. The problem of verifying the validity 
of the above picture is still unsolved in general, and in the situation where the system 
is located in three dimensional space, which is most interesting from a physical point of 
view, there are very few results. However some progress has been made in lower dimensions 
both for finite systems (A^ < oo) and in the thermodynamic limit (A^ t oo). The case of 
one space dimension seems to be fairly well understood [3,4,6,7,20], but see [36]. For 
the case of a two-dimensional space our understanding is more restricted [8,9,10] and in 
higher dimensions only very few results are known [5,24,]. In one and two dimensions 
also models with quasicrystalline behaviour at zero temperature [11,12] have been studied. 
Quite recently a quantum mechanical model showing crystallization has been discussed 
[13,14]. For a very readable survey of current insight into the crystal problem we refer the 
reader to [15]. 

The present contribution is devoted to the analysis of the crystal problem in three di- 
mensional space. We select a special class of interactions for which the analysis of the 
thermodynamical limit of the free energy density describing fluctuations of a system at 
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low temperature around the hypothetical crystalline configuration, which minimizes the 
relative potential energy, is possible by relatively simple and standard methods. As we 
will see, the condition for the minimizing potential energy configurations can be expressed 
mathematically as a condition of positive definiteness of the Hessian matrix for the cor- 
responding potential function. This enables us to rewrite the higher order contributions 
of the Taylor expansion around a crystalline, energy minimizing configuration in terms 
of certain Gaussian integrals. In this picture, non Gaussian fiuctuations around the crys- 
talline structure can be described by a highly nonlinear interacting system of dipoles in 
which the potential describing the interaction between the dipoles in the dipole gas frame- 
work decays rather slowly, but the interaction energy is intcgrablc. A similarity with the 
corresponding problems of lattice systems may be seen from this picture. 

The continuous and the lattice dipole systems with long range like interactions have been 
studied intensively only recently and a variety of deep results describing their high and 
low temperature properties have been derived rigorously. 

We expect that some of the results obtained in [16,17,18,19,20] can be applied to the 
dipole systems which we put in evidence in our study of the formation of crystals in the 
three-dimensional space. 

Restricting further the class of admissible potentials we consider a class of potentials which 
leads to bounded Taylor remainders for the corresponding expansions around crystalline 
ground state. For such potentials the associated dipole system describing fluctuations can 
be interpreted as a grand canonical ensemble gas of charged dipoles living on a Bravais lat- 
tice floo and in thermal equilibrium at temperature where T,. is the temperature of the 
real system. Thus the low temperature behaviour of a real system is described by the high 
temperature behaviour of the corresponding charged dipole system. The activity of the 
associated dipole system is given by — where k is the Boltzmann constant and thus the 
low temperature behaviour of the real system is described by the high temperature/high 
density limit of the corresponding dipole system. Section 3 of the present paper is devoted 
to the study of the thermodynamic limit for the corresponding dipole systems using some 
basic tools of statistical mechanics, e.g. linked cluster expansion and Kirkwood-Salsburg 
identities [22,23]. Although most of the results established there are valid for an ensemble 
of dipoles of arbitrary sizes we shall restrict ourselves to the study of restricted ensemble 
consisting of dipoles of a uniformly bounded size. This cutoff is necessary for obtaining a 
finite contribution to the free energy density of the system we consider. The restriction to 
an ensemble of uniformly bounded dipole is equivalent to the restriction of an admissible 
size of possible deviations of real particles configurations from cristalline configurations. 
Prom the Gibbs measure point of view this means that we have restricted ourselves to the 
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contribution to the free energy density coming from locally finite configurations with uni- 
formly bounded (at large distances) deviations from the (hypothetical) crystalline ground 
state. The problem of large deviations which we called large fluctuation problem is the 
main topic of a forthcoming, second part of this work [34] . An additional discussion of the 
meaning of the restriction to the bounded dipole case is included in Section 4, where also 
a discussion of the compatibility of our assumptions on the interaction potential paper can 
be found. 
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Chapter 2 



Expansion around a crystalline 
ground state. The finite volume case. 

Let us consider a system composed of n classical particles with masses equal to 1 (for sim- 
plicity) and enclosed in the bounded region Aat = {x & \ < N for i = 1,2, 3}. The 
phase space coordinates of the system will be denoted by (p, a;)„ = {{pi,xi), ...,{pn,Xn)), Pi G 
R^, x-i E A^v- The particles interact through two-body potential forces described by the 
potentials V . The class of potential for which our analysis applies will be determined by 
imposing certain conditions which will be listed below. A preliminary discussion of the 
compatibility of these hypotheses is given in section 4 of the present paper. 

Hypothesis 

The potential V is central, i.e. there exists a real- valued function ^ on R\ such that 
V{x) = $(1 X I). The function $ has compact support and is of class at least C ^. The 
potential V is stable in the thermo dynamical sense, i.e., 

3B > 0, Vn : S{{x)n) = ^ V{xi - Xj) > -nB. (2.1) 

l<i<j<n 

The canonical (Maxwell-Boltzmann ensemble) partition function of the system at (inverse) 
temperature /3 > reads: 

ZnW)= [ , f[dpidxiexp[-p£{{p,x)r,)], (2.2) 

where the energy function is 

^((p,x-)„) = E|+ E V{x.-x,). (2.3) 
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Integrating out the momenta we obtain: 

= i^f"^' L ®Udx^ exp[-/3^((x)„)]. (2.4) 

It is natural to try to apply the saddle point method to the study of the low temperature 
{(3 t oo) behaviour of the canonical partition function Z^{i3). The classical Laplace 
theorem tells us that if the minima of the function S{{x)n) are well separated then the 
following asymptotic formula is valid: 

r " 1 

zUP)moo- E e-^''^''"'''^ l^^n^M-PT^iM-xrn (2.5) 

ymin N i=l ^ 



•exp[-/?7^(Xf■^;(x)„)] ^^=1 ^a;,. 

where we have denoted by X™™ = (x^, x^) the minimizing configuration for the energy 
function and the contribution coming from the first factor i.e. from e"^^^'^" dominates 
as /3 t cxD. 

The zero temperature crystal problem in this picture amounts to answering the ques- 
tion whether the energy function E has a minimum at an (almost) regular configuration 
X™™ = (x™, X™). In the finite volume situation we do not expect that the minimiz- 
ing configuration X™™ forms a completely regular structure, because of surface effects. 
To avoid such problems we therefore pass to the thermodynamic limit, i.e., we consider 
the limit X t oo, keeping the density p = n/{2Nf' fixed. It is expected that there 
exists a certain regime of values and p in which the limiting minimizing configurations 
X™™ = (x™, ...,x^) are quite regular globally. A particularly interesting situation is the 
one where the (hypothetical) energy density minimizing configuration X™™ forms a Bravais 
lattice in 3. [| 

Let us recall that a Bravais lattice VL^o in is defined by three linearly independent 
vectors Oi, d-2, % G such that 

^To deal with the case = oo it is convenient to use the energy density function 

ep((a;)oo) = hm g((x)„) ^ 
p ~ fixed 

provided it exists. 
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Let us set = A^nQoo- Assuming n =\Qpf\ wc can, by suitably relabeling and changing 
the variables, rewrite the energy function of n-particles in the following way: 

^ii^)n)^l E V{X-^i + y,-y,) (2.6) 

A, G 

where the variables yx measure the deviation of the A particle from the lattice point 
A e Qat- This can be further represented as: 

£{{x)n) = ^^'■(l^iv) + S'^^^in^ I {yx}), (2.7) 

where 

£-(n^) = ^ E v{x-t,) (2.8) 

A ^ /i 
A, // e 

is the energy of the crystalline configuration fl^ and 

£'''\nN\{yx})^S((x)n)-£'^(nj,) 

= \ E [V{X-iJi-ryx-y,)-V{X-i,)] (2.9) 
A, /X e Jljv 
A^/i 

measures the energy deviation of the real configuration {x)n from the energy of the crys- 
talline configuration Ojv- 

The crystalline energy E'^'^ depends only on the parameters ai, a2, spanning our Bravais 
lattice. Therefore one can easily find an equation on ai, 0.2 03 that yields a stationary point 
for the crystalline energy £'^'^{^1^) (see e.g., [21]). As an example, in the case of central 
potential forces described by the function $ as in Hypothesis the corresponding equation 
reads: 

E lY'i = 0, ^,J = 1,2,3, (2.10) 



n 7^ 



which is a set of 9 equations with 9 unknowns. This equation, together with the condition 
that the relevant Hessian is positive, then gives a local minimum for £^"^''(Q„). 



6 



We find that for technical reasons it is more convenient to work with a periodic version 
of the given finite Bravais lattice, henceforth we shall understand under fl^ this periodic 
version. 

The following hypothesis plays a crucial role in the present attempt to understand the for- 
mation of crystals in the low-temperature/high density region in three dimensional space. 



Hypothesis 1 

(Hla). There exists a number pcr > such that for all p > per there exists a Bravais 
lattice solutions f^Ar of the minimizing problem for the corresponding £'^{0,^). Moreover 
this holds uniformly in N at least for sufficiently large A^. 

(Hl/3). The deviation energy function £'^^'"{yL]s[ \ d/njAea^) has a stationary point at 
{y\ = 0, A e Ojv} which is a local minimum for sufficiently large values of A^. 

Part /3 of Hypothesis 1 can be derived from the condition. 

< (y, [DV]t/) ^ ^ E {y^-y^^[D'V]{\-p){y^-y,)) (2.11) 

A,^er2jv 

where [L'^l^](A — p) is the Hessian of V evaluated at the lattice point (A — /i), see e.g., 
[21]. 

Condition (2.11) can be used to deduce the existence of a certain functional integral rep- 
resentation of the quantities of interest. For this reason we define the matrix: 

A.(A,,) (2.12) 

Then we have the following equality: 

I E {yx-y,:[D'V]{X-p){y,-y,))= ^ (yx, Aj,{X, p)y,) (2.13) 

From the hypothesis HI/? it follows that the matrix An is strictly positive definite on the 
space R^l^^l (with \^In \ the number of points in ^at). On the Borel cr-algebra of sets 
of i?3|f^iv| Qg^j^ therefore define the Gaussian measure p% with mean equal to zero and 
covariance matrix Cn{X,p) equal to the inverse of the matrix An{X,p), i.e., 

/ijain.i VAr(2/i,...,y|n.|) exp[i ^ {a,, y,)] ^ E%{e'^^'^) (2.14) 



= exp[-^(Q;, A^^a)] = exp[-^(Q;, Cno)] 
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where 



a = (cti, ...,a\nj^\), (2.15) 



{a,y)= («a,|/a), (2.16) 

{a,CNa)= ("A, J2 Cn{\iJ')oii^), (2.17) 

and similarly for (a, 74^^(5). From the periodic lattice Qn and the infinitely extended lattice 
Jloo we form the corresponding dual lattices Fjv and Too respectively. The corresponding 
dual groups or Brillouin zones are denoted by Qjv and ft^ respectively and are defined by 
= [-N, TVJVrjv, doc = -R Vroo- We have: 

I E (^A - Vi., [D^V]{X - i^){yx - y^)) = ^ ^^(i^(0) - iAr(p))yiv(p)(2.18) 
A, At e Qjv ^^^^ 



where 



and 



A^ip)^ Y: [D'V]{X)e^^' (2.19) 

Aef2jv,A^0 



yr,{p) ^ Y 2/Ae^"'- (2-20) 

AeQjv 



For the infinitely extended lattice fl^o the corresponding formulae are: 



1 
2 



^ (y^ - y^, [D'V]{X - f,){y, - y,)) = / y(p)(ioo(0) - AM)y{p)dp{2. 21) 



A^/i 

where 



= E ?/Ae^^^ (2.22) 
AeHoo 



ioo(p)= E [D'V]{\)e^^^ (2.23) 

A G fioo 

Ay^O 
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We remark that, due to the compact support of D'^V{X) and the fact that fl^ is compact, 
Aoo{p) is well-defined, and is the pointwise limit of A^{p) as N ^ oo. Similarly, the 
right hand side of (2.22) is well-defined for all y G ^2(^00)- Hence the equalities (2.21), 
established first for a dense subset of L2{floo), extend to the whole of L2(0oo)- 

Let us define the functional Sq^ on the dense subset of L2{floo) consisting of sequences 
{yx}, A e Qoo with compact support, by the formula (2.9) with N — 00, i.e., by 

Sf\n^\{y,},^^J^^ Y: [V{X-^i + y,-y,)-V{X-^i)]. (2.24) 

A, /i G r^oo 

Then we have, using formula (2.20), the following: 

Lemma 2.1 Assume that the point {y\ = 0}, A G Qoo is a stationary point of Sq^- 
Then it is a local minimum for Sq^ if Aoo{0) — Aoo{p) is a positive definite matrix for all 
p G f^oo- The corresponding statements hold for S^^^ i.e., if {yx}, A G is a stationary 
point for £[1^^ then it is a local minimum if Am{Q) — ^jv(p) is a positive definite matrix for 
all p & Clf^. 

By Fourier analysis, similarly as above, we find the following representation for the inverse 
matrices {An{\, fJ'))~^ and Aoo{\,fi)- 

Cn{X -pi) = (^iv)(A,/x) = e*"^^"'^n^iv(0) - iiv(p))-' (2.25) 

and 

C^i\-fi) = {A^'){\,lx)= f e^f(^-'^)(ioo(0)-ioo(p))-'c?p. (2.26) 
The following hypothesis will play an essential role in the following analysis. 

Hypothesis 2 

(B.2a). There exists a constant c, independent of N, such that A^ip) — An{0) > cp^ -\- 

^Nip^) as p ^ 0. 

(H2/3). i^(p) = i^(0) iffp-0. 

Hypothesis 2 ensures that the functions Cjv(A — //) and Coo (A ~ A*) are well-defined and 
moreover Coo (A — A*) decays at least as p^^^ for | A — /i |— > 00 . The set of potentials for 
which (H2q;) and (H2/3) are satisfied is surprisingly rich, we refer to the work [21] and to 
section 4 of the present paper. 
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The previously introduced Gaussian measure dii% has the following density with respect 
to the Lebesgue measure 



N 



det[(^^iv(A,/x))A,^en^] ^exp[-^ ^ {yx, An{X, |^)y^^)] (2.27) 



Expanding around the minimizing configuration 0,^, using the Gaussian measure integra- 
tion and making a simple change of variables, we can represent the canonical partition 
function in the following way: 

Z%{/3) = (27r)3"/2^-='" ■ Z'^in^) ■ ZN ■ Z-fj{(3), (2.28) 

where we have introduced 

• the crystalline partition function: 

Z^'Xn^) = exp[-[pS"\nN)] (2.29) 

• the fluctuation partition function: 



ZjjiP) ^ E%{x{{\ yx \< N}xen^)) exp[-/37e(n^ | {^}xen.)] (2-30) 

where 



niQN I {yx}xen^) = S{{x)n (2.31) 
-I E {yx-y„D'V{X-i,){yx-y,))-S'^^{nr,) 

is the Taylor remainder of the Taylor expansion of the energy S{{x)n) of the real 
configuration (a;)„ around the local minimum given by fl^ 

• the normalisation factor for the measure 

ZN = det(;^(Ajv(A - fi)x,f^ea^y^^ (2.32) 

ZTT 

The free energy density of a system enclosed within the region Ajv (with periodic boundary 
conditions) is defined by: 

pp^{P) = ^\nZ^{P) (2.33) 
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Substituting (2.28) into (2.33) we obtain: 

Q 1 r) 

PPn{P) = -\^P - 3pln/5 + (iP^^i^N I (5) + (5~ZM + (5pn{(5) (2.34) 

where 

Pn{^n I P) = -^-r^lnexp[-/3£^'-(f)^r)] = -^^T (^.35) 

p I Ajv I I Ajv I 

is the crystaUine free energy density; 

= T^^r^lndet((A^%^A,/x e Qjv) (2.36) 
is a numerical constant and 

^^^^^ ^ i^"^^""^^^^^"^' ^}Aen.)exp[-/57e(0^ | j^j)]] (2.37) 
is the free energy density of nongaussian fluctuations around Qtv- 

A simple Fourier analysis is used to control the thermodynamic limit of the quantity zjq. 
We have 

~ZN = 4o^lndet((A^%^A,A^ e Q^) (2.38) 

zp I Ajv I ZTT 

-^arr-\ ^ in(iiv(o)-iiv(p)). 

2^ I I pen. 

As A?" — > oo this tends to: 

-^oo = / dp tr ln[ioo(0) - ioo(p)], (2.39) 

2/3 I A I 

which is finite, e.g., under the assumption (H2q;). 

The most difficult element of our analysis concerns the problem of the existence of the 
thermodynamic limit of the free energy density of fluctuations around a crystalline ground 
state. Taking into account the formula 

7e(l^^ I {^},enj = E 7^^.(A-/.|^^), (2.40) 



where 



7^^,(A - ^Ji I y^-^) ^ V{\ + - - /S-'/'z/m) (2-41) 

_V(A - ;,) - ^^^^ PV](A - /.)(^^)), 
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we see that in the field theoretical picture given by the Gaussian integral representation 
(2.30) the partition function Z-y can be viewed as the partition function of a gas of dipoles 
that sit on the lattice and interact by a nonpolynomial and in general nonlocal po- 
tential given by (2.40). With some additional simplifying assumptions the problem of 
the thermodynamic limit of the free energy density and also for the corresponding Gibbs 
states describing the associated dipolc systems will be studied in the next section by an 
application of techniques of statistical mechanics e.g. linked cluster expansion and analysis 
of the Kirkwood-Salzburg identities. 
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Chapter 3 



Infinite volume limit. Bounded 



dipole length case. 



In this section we shall discuss the thermodynamic limits of the free energy density and 
the Gibbs states of the associated dipole systems describing bounded fluctuations of the 
real configurations around the crystalline ground state given by a Bravais lattice Qoo- 

3.1 A convergent expansion for the dipole partition function In order to simplify 
the analysis as much as possible the following additional hypothesis will be introduced. 

Hypothesis 3. There exists a family {dfix)x(zQ^ of complex measures on such that the 
following representation for the corresponding Taylor remainders defined by (2.42) hold 
for all A — /i: 



In addition we assume that dfi\{a) = dfi\{—a) so that TZn^i^ — /i | x) is real. We also 
make the following additional assumptions on the measure appearing in (3-1)^. 

(H3)a). Jd\fix \ (tt) I a |"< oo, for all A and all n G A?" [] 

(H3/3). Jd\fix \ (a)e^" , for some ^ > 0, and all A e floo. A has to be sufficiently large 
(see below). 

(H37). For all A, supp dfix is compact in R^. 

The following trivial chain of implications holds: (H37) =^ {H3/3) =^ {H3a). Which one 
of the hypothesis a, (3 and 7 is needed depends on the particular situation we investigate 
(see below). 

As a preliminary step in our analysis we remove the characteristic function restricting the 
^In fact we need only the existence of the first moments (see Proposition 3-2 below) 




(3.1) 
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possible values of the size of fluctuations {yx} in (2.30) defining 

Zff = E%{exp[-pn{nN I {p-'^'yx}xenj). (3.2) 

Using the formula 

E%{l[{e'^~''""^''^ ■ e-^^-^'^'^^y-^) = exp[-i ^ V,v(A./..a. | A,/x,-a,-)] (3-3) 

where 

VNiKl^iOii I Xjl^jaj) = ai{CN{K - Aj) + CNil^i - j^j) (3.4) 

-CAf(Ai - /Xj) - CAf(/Xj - Ai))Q;j 

and taking into account H4 and H6 we obtain the following expansion of Zf^: 

^^ = E^ E /®r=irf/^A,-«(«o (3.5) 

I Aj — //j I < i? 
^^P[~^ E ^iv(aiAi/Xj I ajXjiij)] . 

P l<i,j<n 

From the very simple inequality 

\ E%(eMiT.o^xyx] I<1 (3.6) 

X 

we deduce the following estimate on Zf^: 

Z% < exp[| L>jv(l) |] (3.7) 



where 



II* = sup { / I /xa |}, 



and I -Djv(l) | is the cardinality of the set Qjv x ^n- Estimate (3.7) shows that the 
expansion (3.6) is absolutely convergent for any /3 > 0. For the dipole free energy density 
p%P) =1 -Div(l) InZ^' we obtain the upper bound: 

pUP) < Pl^* (3.8) 

uniformly in N. 

Although some of the results below hold for a dipole systems described by the dipole 
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partition function (3-3) it follows from the estimate (3.8) that in order to ensure the 
finiteness of the limit 

lim — — rlnZir 

(see for (2.37)): without further assumptions on the interaction potential V we have to 
restrict the admissible length of dipoles. Let it! be some fixed positive number. Then we 
define a bounded length dipole system by the following partition function 

Z'M = E%{e^p -PKi^N I {p-hx}) (3.9) 
where the restricted Taylor remainder R'^ is given by: 

Ki^nlirhx)^ E Rn{x - ^i \ ^^^^^) (3.10) 

I Aj — I < 

From now on an allowed maximal size of dipole given by a number R will be fixed. The 
following sort of stability also comes from estimate (3.7) in a simple way: 

n 

E Viv(Q;jAiA*i I XjlJ'jOij) > -2E«i(C'iv(0) - C'iv(Ai - iJ,i))ai. (3.11) 

l<i^j<n i=l 

The expansion like (3.6) for has the interpretation as a grand canonical partition 
function of a system of dipoles of length bounded by R on the lattice f^^v in the thermal 
equilibrium at inverse temperature l/f3. Because of the crucial stability estimate (3.11) 
the standard tools of classical statistical mechanics, such as cluster expansions, Kirkwood- 
Salzburg type analysis [22] can be used to study the corresponding thermodynamical limit. 
In order to simplify our notation the following abbreviations will be used: 

D§{K) = D§{1) X ... X D§{K) (3.12) 
= {(Ai,//i,q;i), (A;^, /XftT, ax) | Xi^pii £ ^n, | Aj - /Xj |< i?, aj G R^} 

= K-dipoles configuration space. 
Furthermore we denote by 

D§{K) = D^{1) X ... X D§{K) (see (3.10)) (3.13) 

the X-dipoles restricted configuration space. 

For a generic u; e D§{K) the following abbreviations will also be employed 

UJ = ((Ai, 111, Q!i), (A„, l^n: "n)) = (A, /X, «)„, -0^(1) 3 (Aj, /Xj, Qlj) = d{i). (3.14) 
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/ 4«- ^ 



J2 j df^Xi-miai) 

K — fj'i I < -R 



(3.15) 



/<(l,...,n)... = /4(rz).../4(l)- 



/ ^^(0- 



E 

Xi — Hi \< R 



OLi 



(3.16) 
(3.17) 



£iv((Qi, A,//)n) = ^ Viv(Qii, Ai,//i I ajAj/Zj) (3.18) 

l<i J<n 

£oo((«, A,/x)n) = E T^oo(aiAj/ij I ajAj/Xj) (3.19) 

l<i,j<n 

£:jv((a,A,/i)„| (q;',A',//OJ (3.20) 
= £^Ar((Q;, A, )„ U (a', A', //')m) - £^Ar((Q;, A, //)„) - £^Ar((Q;', A', //')m) 

and similarly for £^oo- 

Hence we can rewrite (3.6) for in the following way 

^^' = E^^ / <(l,...,n)exp[-l^^(rf(l),...,c^(n))] (3.21) 

n>0 J P 

3.2 The linked cluster expansion for the free energy density. 

Using a well known trick of Mayer [22,23] we can derive the linked cluster expansion for 
the free energy density for a finite volume system. For this purpose we write 

exp[--£^,(d(l), d{n)) = n ^A^(^) 11 exp[-- V^(d(i) | d{j))] (3.22) 

P 1=1 l<i^j<n P 

-\{VN{^) n [(exp-V^(ci(z)|ci(j))]-l) + l] 

n 

= ^^iv(^) E M{d{ieT)) 
j=i rc{i,...,n} 

where 

M{d{i e r)) ^ n {^M-\vN{d{i) I d{j)] - 1) (3.23) 
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and 

V^it) = exp[-^a^{CN{0) - CM{X^ - (3.24) 
is the corresponding vertex function. 

We denote by the set of all connected linear n-graphs that can be built on the set 
D^{n). Putting the expansion (3.23) into (3.22) and resumming we obtain the linked 
cluster expansion in the following standard form: 

Z^^^expE (-/?)" -6^^] (3.25) 

n>l 

where 

b%= i: &^(r) (3.26) 

and the contribution b^{T) from the graph T G is given by 

^'^^(r) = ^ / 4(1)... / 4(n) n [exp-^V^(0 - 1] n (3.27) 

where we have denoted the set of lines of a given F e by C{V), by V{V) the set of 
vertices and taking Z e F we denote by Iq the initial and resp. by 1^ the endpoint of I and 
then we define V nQ) = N{h \ U)- 

Prom the positive definiteness of Cat it follows that Cn{^) — Cn{ij) > for any ji e 
Jljv and uniformly in N. Therefore we can estimate the vertex function contribution to 
(3.26) by: | Y{i^v{T)^N{i) |< L Our first result for the dipole free energy density in the 
thermodynamic limit A'^ ^ oo is the following. 

Proposition 3.1 Assume that the hypotheses H2, H3 and H3j are all valid. Then: 

yn^T e gjh^moo \ b'^NiT) \< oo for any /3eR+-{0}. (3.28) 

Proof: 

Any connected, linear n-graph T G Qn contains at least one spanning tree st(r). Let 
I C{r) 1= k + s, where A; = n — 1 is the number of lines forming st(r) and s =\ C{r) \ — n + 1 
is the number of lines in r\st(r). The contribution coming from the line / ^ st{r) can be 
estimated by a "sup argument": 

sup I exp[iv-jv(0] - 1 |< | ai, - ai^ \ exp\j3-Hai,ai^C *] (3.29) 

I P 
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where 



Setting 

it follows 

where 

and 

b^(\ str \) = 



C* = sup(sup II CM II) 

N M 



a* = sup sup{| a II q; e supp dfix} 



b^{r)\<E{l3)X{\stT 



E = AI3-^C*a*^ exp[/3-Ma*2C*] 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



Proceeding further we obtain 

1 



j 4 I (l,...,n) I n I (exp[--y(ci(0 I d(^ + 1))] - 1) I . (3.34) 



C(l stV I) |< -jr^E'^"-'^ / I 4(1,-,^) I n I VN{d{i) I d(i + l) I (3.35) 

I ^N\^) I i=l 



where now 



£;' = exp[-a*2<^*]. 



(3.36) 



Let Tp be a translation by p G VL^. Observing that the equality TpVjv(AyU« | X'/j^'a') = 
Vn{X + p, 11 + p,a I X' + p, p' + p,a') = VN{Xpa \ Xpa) holds, we can estimate 

/n— 1 
|4(l,...,n)| n \VN{d{i)\d{i + l))\ (3.37) 
1=1 



< o(/3) E E 



E 

An, /^n 



A2 — A*2 I < -R I A„ — //„ I < it! 



I /"Ai-Ati I (ttl) ■ / I /^A2-M2 I ("2)--- / d I 

/-'■Ti /-'■n I iS^n) 



Voo(0, Ai,ai I A2Ai2a2)---l^oo(A 



<oy) E 



AiSOfl 



E {j d \ f^Xil I /^A2-/.2 I («2) 



A2 - A*2 |< -R 

A2, //2 G ^^oo 



^"00(0, Ai, «! I 0:2, //2, 0:2)} 



n-1 
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where 0{^) and 0'{^) are numerical constants depending possibly of R but not of A^. We 
finish the proof if we show 

E E j d \ I^Xil Md I /^A2-M2 I M I (3.38) 

^ A2, /i2 £ ^^oo 
I A2 - //2 I < -R- 

•Voo(0, Ai,q;i I A2,A*2,Q;2) |< 00 

This amounts to study the decay properties of \4o(0, Ai | X2, 1^2) = C'oo(A2) — C'oo(/i2) + 
Coo(Ai - fi2) - C'oo(Ai - A2) as I A2 IT 00. 
Let us define 

C,{X)^C^{X)-C^{X-pi). (3.39) 
Taking the Fourier transform of (3.40), using (2.21) and that Coo(A) = Coo(— A) we obtain 

= (1 - ^cosp//)(Aoo(0) - AM)''- (3-40) 

Prom the hypothesis (H2/3) it follows that Cij,{p) stays bounded as p — > whenever // 7^ 0. 
Therefore by the Fourier- Tauberian theorem it follows that C^(A) has to decay at least as 
I A as I A II cxo. Therefore taking n' ^ 0, we conclude that the function C'oo,ft(A) = 
^n{X) — C^{X — n') has to decay at least as fast as | A 1^"^ as | A If 00 which is intcgrable 
over R^. Let us apply this information to the problem (3.39). Using the definitions of 
Coo and Coo,^ 

E E I ^oo(0//iQ;i I X2H2a2) \ (3.41) 

/i2, A2 t iioo 
I A2 - //2 I < -R 

<l «! I ■ I "2 I "E E I C'oo,/.i(A2) - C'oo,^i(/i2) I 

/X2, A2 G r^oo 

I A2 - /i2 |< -R 
<| «1 I • I «2 I •0"(1) • I ^oo,m(A2) |< 00 

for some A* 7^ 0. | | 

Remark 

The above Prop. 3.1 can also be proven, essentially by the same method, replacing as- 
sumption H37 by H3/9. 

Let us define the virial coefficients fe^(r) and the limit t 00 as the limit for f 00 
of 6jv(r) as given by (3.28). Moreover by the translational invariance we impose the 
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restriction on F that at least one of it vertex has to be located at some (0, A) where 
A G Or- That the limit limTvtoo ^7v(r) exists follows from the simple observation that for 
a fixed F, the sequence (6^(F))jv is a Cauchy sequence in N. 

Let us remark that the limiting virial coefficient 



X n Vi^ (3-42) 

iev(r) 

are analityc functions in for /3 7^ and moreover lim^_»oo &^(F) = 0. Prom the Taylor 
expansion it follows that: 



^ ^ ^\c(r)\ , , ^\v(r)\ , (3.43) 



ki,...,k\c{T)\ r-j,. . . ,r\v{T)\ 
ki>0 n > 



\c{r)\ \vir)\ 



Ail e Or 
(0,/ii)e£(F) 



• n [exp{~vM - - o{k))] n ^(^) 

/e/:(r) ^ iey(r) 

from which it follows that also 

hm 'T, J ^ (3.44) 

Since the total number of n-graphs F restricted as above is finite, we conclude that the 
virial coefficients 6^ are analytic functions of for /3 7^ and moreover limjS^oob^oiP) = 
and also 

lini = 



20 



for any value of M. 

A very convenient way to describe the virial coefficients in the hmit N ^ oo has been given 
by Brydges and Federbush [24] . We adopt their method to the case of the dipole systems to 
give the hnked cluster expansion for Poo{P)- Towards this goal let us denote by Vn the set 
of tree functions 77, defined on the set {1, ...,n}, i.e., 77 e T„ 77 : {1, — > {1, ...,n} 
and r]{i) < i for any i. Denote by = (si, s^-i) G [0,1]®""^ a sequence of 

interpolating arguments and define corresponding interpolating energies by the following 
inductive process 

E oo((A, //,«)„) = £^oo((A,Ai,a)n), (3.45) 

^ oo,i((A, n, a)n) = (1 - Si)S^{{X, n, a)n \ (A,, i^i, a,)) + Si£oo,i-i((A, yu, a)n} 

S oo,n-l((A, /i, a)n) = S^{{s)n-l). 

Define a function f{r), (s)„_i) by 

fiv, (s)n-i) = Ili=i Si-iSi-2-Sr,ii), for n > 2 
f{v,si) =1 

where Sj_i...s„(j) is 1 if ?7(i) = 1. The n-th virial coefficient can be written in the 
following compact form 

= E / di^)n~i J d^iX,fi,aU (3.47) 



n— 1 -1 
fiV, {S)n-l)) n + 1)) I divm eM-l^£Zi{s)n-l)] 



i=l 

where 



/9 



j d^{\,fx,a)n... = E y"c^/^w("i) E j dixx2-tMi ■ ■ ■ (3-48) 

A2, ^2 £ ^00 
I A2 - /i2 I < -R 



The following estimate is well known (see [24]) 

Y: [ d{s)r,.J{r), (s)n_i) < e"-^ (3.49) 

The induction steps defined in (3.46) involve the convex sums of energies only, therefore 
the positive definiteness of the intermediate and final energy function is preserved in this 
process. This observation and estimate (3.50) lead to the following bound on the n-th 
virial coefficient 

1 p"~^ 

I c i< Y^i^ II iir^ (3.50) 
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where 

II Voo ||i=| Or\ sup{ d\ nx\ {a) \ a |) (3.51) 

■ sup ( X] / ^ I A^A'-M' I (q^') I Q^' 1 1 ^00(0, A I A', //') I 

A',//' e Qoo 
I A' -//' |< 

where 

Voo(A, /X I A', = Coo(A - A') + Coo(/^ - 1^') - Coo(A - - Coo(/x - A'). (3.52) 
In this way we have proved: 

Proposition 3.2 Assume hypotheses H2, H3 and H3j are satisfied and moreover that 
II ^00 ||i< then the virial expansion 

PM = ■T.-T. I d(s)n-i [ d^(X, I,, a)° (3.53) 

n>l ^ r?eT„ 
n— 1 -1 

/(r/, (5)„_0 ■ n ^oo{d{z + 1) I d{l^{^)) exp[--SZ{{s)r,-i)] 
i=i P 

is absolutely convergent for any /3 e C such that Re(5 > 0. 
RemEirk 

Define the Borel transform Vp{Z) generating function for the virial coefficients by 

. E '-^>°7-'^' (3.54) 

n>l 

Prom the estimate (3.50) it follows that VpiQ is an entire analytic function on C in 

for any j3 G -R'''\{0}. Therefore applying the inverse Borel transform to Vp we obtain a 
function Vpi^z) = {B~^V^ foi' which an expansion in powers of z coincides with the 

virial expansion (3.54) for 2; = 1 and moreover the expansion is a Borel summable as /? J, 0: 

W\.=i = E + Oiv(/?^+^) (3.55) 

n>l 

with I Ojv(/?'^+^) \<\Cn\P for some constant Cn > 0. 
Remcirk 

The Borel summability of the high-temperature expansion in the classical statistical me- 
chanics has been studied in [26]. To prove Borel summability of the corresponding low 
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temperature expansion (3.54) one needs to expand the corresponding Mayer kernels in 
powers of and then reformulate the whole expansion (3.54) in powers of 

3.3 The associated infinite volume Gibbs states. 

In this section we shall construct the infinite volume Gibbs states describing dipole systems 
fulfiUing all hypotheses HO - H3 above. It is well known that the corresponding Gibbs states 
are determined by their correlation functions, see e.g. [22]. Therefore we concentrate on 
them in the following. 

The finite volume correlation functions are given by: 

n 

pjv((A, /.,«)„) = {-f3rE^il[exp[tf3-'/'aj{y,.-y^.)] (3.56) 
where the expectation denoted as is defined as 

^ E%{eM-m^N\{(5-'''yx]xen. ^ ^ 

To be more general we consider also the possible influence of the external dipole config- 
urations on the thermodynamic limits of the corresponding Gibbs states. Towards this 
goal, let us consider a generic point u e (Qoo x ^oo x -R^)°° = X>p(oo), where I>p(oo) is 
defined as a configuration of dipoles such that to every point (A, //) e Qqo x ^oo there is 
inserted at most p dipoles with p is some integer. It follows from the finite length dipole 
approximation that for every oj G Vp[oo) there is only a finite number of dipoles that can 
share the same end point A e Q^oo and which contribute to the total energy in a nontrivial 
way. 

The restriction of agiven dipole configuration a; e ^p{oo) to the set (QooXJlooX-R^)\{JlArX 
Qjv X il} = will be denoted by uj{N'^). The associated, finite volume, conditional 
Gibbs states are described by their conditional correlation functions p*^ which are defined 
for uj e Jlp(oo) by 

p5^((a,A,/.)„) = (-/3)".£;^(n^^^exp[^r^«,(i/,^. -y^J] 
•exp-/37^(^]^|■) -exp -(3J'm{^{N^))), 



(3.58) 



where 



J^,,{u{Nc))^ E 7e(A'-/.'|^^^) (3.59) 

gives the energy of a given configuration of dipoles (A, /i,Q;)„ inside Ajv with an external 
configuration uj{N^) realizing a particular boundary condition. Due to the decay properties 
of Voo and the assumption uj e X'p(oo), the energy in (3.59) is finite. 
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The conditional partition function used in (3.59) is given by the formula: 

^;°^(exp-/3^^(a;(7V^))) ^'-'"^ 
Performing the following, (complex) shift transformation 

in formula (3.59) and respectively 

yx-Vn^ y\- iC^iyx - Vm) -Vh- iCNil^ - l^i) (3.62) 

in the formula (3.57) and calculating the corresponding Radon-Nikodym derivatives we 
obtain the following identities 

p%{{a, A, /i)0 = {-f3rexp[-f3-'SMm\ f^, «)n U u;{N^))] 

•^^^(n^=2e*''^"'^'''^"''''^^exp-/37e(f7iv|e-^"'^-(^WI-) _ i)) 

exp -pj^Ni^{N^)\e-f^~'^-m^ - 1), 

(3.63) 

where 

^Ar(a;(Arc)|e-^-^^-W)l-) - 1) = J2 / ^i7(«'A_^)e"^-''^^\e-^"^-('^(^)IK-..^.^) _ i) 

\X-IJ,\<R 

and 

Pjv((q;, A,/i)n) = exp -p-^SN(d(l)\(a, X, ii)n) 

eUU]=2 e'^'^'^^^'^^ry^i^ exp -/37^(Qiv|(e-^"'^(<^)l•) - 1)) (3.64) 

where 

7^(^]^,|e-^-^^-('^«l•) - 1) = E / diix-Mx-,) ■ e^"-'^('^\e-^"^^-('^(^)IK-..^'.M) _ i) 

\\-ii\<R 

in which, after an expansion, we recognize the well known Kirkwood-Salzburg identities, 
see e.g. [22]. Using the method of a dual pair of Banach spaces as explained in [25], a 
large region of admissible (3 will be determined in which the rigorous comparison analysis 
of (3.64) with (3.65) is possible. 

For any ^ > introduce a Banach space E^^ consisting of all sequences of functions 
(/„((A,//, q;)„)„=i,2,...) such that for any n, any A,// e Qoo the functions f\,^j,{{a)n) = 



24 



/„((A, «)„) are measurable and such that the norm || (/„)„ ||^= sup„ ^^"ess sup(-;^^^^^)^ | 
fn{{X, ^,a)n) I is finite. Several bounded linear operators will be defined to act in the 
space E^. The infinite volume Kirkwood-Salsburg operator K^{(3) is defined by 

(^oo((/m)m))n((A, /i, a)„) = exp[-^£:oo((Ai, /Xi, Oi)] | (A, n, a)n) (3.65) 

If mi 

• 5: — / d^(A',//',a')mn(exp[-^^oo((Ai,//i,ai) | {K,fi[,a[))] - 1] 

m>0 i=l ^ 

/n+^_l((A,^, (A', fl',a')m) 

where 

(A, /X, Q;)*_i = ((A2, /X2, Q!2), (A„, i^n, an))- (3.66) 

The finite volume, unconditional Kirkwood-Salsburg operator is defined by: 

(i^iv((/m)m))n((A, /I, «)„) = exp[--^£iv((Ai, ^i, tti) | (A, /X, «)„)] (3.67) 

^ — / 4(A',/.',a')mn(^^P[-75^iv((Ai,/.i,aO | (A^, /x^, a^))] - 1) 
m\ J r 1 P 



m>0 ""■ " i=l 

\l ..I _/ 



/„+^_l((A, «)*_!, (A', //', Q;')m)- 



The energy factor operators are defined by: 

) = Xiv((A, //)„) exp[--£:oo((Ai, /Xi, a,) | a;(7V-))] (3.68) 

■{En ■ f)n{{Kl^,a)n) 

where 

N N I 1 (Aj, itj) e Qat X ^AT for all i 

Xiv(A,/.)„)^ ^ ^ (3.69) 
otherwise . 



E^{{fm)m)n{{\ 1^, a)n) = Xiv((A, /x)n) exp[-^£7v((Ai, /ii, fti) | (A, fi, «)„)] (3.70) 

/n( (A, //,«)„) 

EN{{fm)m)n{{\ a)n) = Xiv((A, //)„) exp[-^£:oo((Ai, //i, Qii) | (A, «)„)] (3.71) 

•/„((A, //,«)„) 

The projections YIn defined by 

n 

^N{{fm)m)n{{K lAn) = H Xn{,{.K l^i))fn{{\ l^)n) (3-72) 

i=l 
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The finite volume conditional Kirkwood-Salsburg operator Kf^ is defined as follows: 

K%^E%-K^- n^. (3.73) 

For further use we also define the operator K^o through the relation 

K^ = En-Koc- (3.74) 

Also the finite volume version of is defined by 

K^ = E^- k^. (3.75) 

Prom the positive definiteness of it follows that we can split the set {0,^0 x ^oo xR^)®'^ 
as U^=iE,", where = {(A, //,«)„ e {Q^ x x R")^^ \ E?^,- V^{\fXiai \ XjfXjaj) > 
—2 Er=i (C'Ar(O) — Cn{K — A*i))}- Let then 77" denote the characteristic function of E" 
and let 9'j — rf- / Ylj^^^rf- . Let 5"^ be defined on functions of n-varibles (A, //,«)„ as the 
circular permutations of /c-steps on the complex arguments (A, /x, a) of these functions. 
Then the index juggling operator J^o is defined as 

n 

Joo{{fm)mU{\f^,a)r,) = ^ -^fc ((^."((A, «)n)/n((A, /X, «) J (3.76) 

k=l 

In a similar way we define the Ruellc index juggling operator J'n which selects a labelling 
of coordinates in such a way that £jv((Ai, /^i, ai) | (A„, «„)) > 0. 

Finally the following vectors are defined 

'^- = (exp[-^£oo((Ai,/.i,ai) I (A,/.,a)iUcu(iV^)),0,...) (3.77) 



aN = (exp[-^£:Ar((Ai,/ii,Q;i) | (Ai, /xi, cti))], 0, ). (3.78) 

With this notation the equalities (3.64) and (3.65) can be rewritten in the following way 

{{pN)n) = (-/3)(•^^^ • Jooi^^ • ^N){{p%)n) + hmNa^ (3.79) 
and respectively 

((Pjv)n) = -Pi^N ■ Jn ■ Kr,IlN){{pN)n) - mNO^N- (3.80) 

We expect that the infinite volume limits ((poo)n) of {{pN)n) should fulfill the following, 
infinite-volume Kirkwood-Salsburg identities: 

((Poo)n) = -PJooKoo{{Poo)n) - (3.81) 
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and in fact we will show this for the infinite volume limits of the conditional correlation 
functions, providing they exist. 

The rigorous, comparison analysis of the identities (3.80), (3.81) and (3.82), based on the 
use of methods of the dual pair of Banach spaces as explained in [25] will be presented 
below. The conclusions obtained in this way seem to be much richer then the contraction 
map principle (see, i.e., [22]) ordinary used for such an analysis. 

It is necessary to introduce the additional (technical) hypothesis for making our analysis 
complete. 

Hypothesis 4 

(H4) There exists a number possibly depending on p and R such that for any ^ > ^* 
the family (p^((A, q;)„)„)jv is weakly -* precompact in the space for any uj e Dp(oo). 

Note that for the case u — (f) and for the correlation functions (3.57) H4 is fulfilled due to 
estimate (3.7) and the Banach- Aloglou theorem. For a general u e -Dp(oo) the necessary 
estimates are still to be proved. 

To formulate the main result of our analysis let us define (7^(/oo-f^oo(/9)) the spectral set 
of the operator J^K^{I3) in the corresponding space E^^ (see below) and let us define the 
set as: 

H^ = {l3eC\Rel3>Q- | /? |< ^ and (-/?)"! ^ (x^(Jooi^oo (/?))} (3.82) 

Theorem 3.3 Let us assume that the hypotheses H0-H3, H3j and H4 are all valid. 

Then there exists a number (possibly depending of p and R) such that for any (3 e H^/^ 
there exists a unique thermodynamic limit p'^ of pff{{\, p,a)n) which does not depend in 
particular on uj ^ Dp[oo) and is equal to the thermodynamic limit p^o of {pn{{\, p,Q)n)) 
which also exists as a unique limit as N ] oo. The limits should he understood as limits 
in the locally uniform, componentwise topology of the space E^. The limiting correlation 
functions Poo depend analytically on fi, provided (3 G H^i , they are translationally invariant 
and posses a cluster decomposition property with respect to the translations offl^ x fl^. 

Before passing to the proof of this theorem let us note some simple consequences which 
seem to be relevant for the discussion of the problem of the convergence of the linked 
cluster expansion worked out in this section. For this goal let us define 

C{P) = sup (I sup{| a I I q; e supp dpx})' I Or \ \ (3.83) 
ijd^ I {X',p',p') I I exp[-i£oo((0,A,a) | (A', /.',«')] - 1 I) 
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and 



C {(5) = sup sup (I sup{| a I I a e supp dux})- \ Or \ (3.84) 
N xeOR 

( / I (A',/.',/.') I I eM-^^NmX,a) I (A', //',«')] - 1 |). 

It follows from our hypothesis that these quantities are finite (see the proof of Prop. 3.1). 
Let us also define 

= {(3 eC\ Re/3 >0 and \ (3 \< min{| ^ \,C-\(3)e-^ < 1}} (3.85) 

and 

Hf = {f3 e C \ Ref3 > and | /3 |< min{| ^ \,C'-\P)e-^ < 1}}. (3.86) 



Corollary 3.4 Let us assume that all the assumptions of Theorem 3-3 are fulfilled. Then 
C H^{j3) and C H^{j3), therefore all conclusions of Theorem 3-3 are true for 
P e i^l (respectively j3 e )■ 

Proof: 

It is very easy to check that for (3 e the corresponding Kirkwood-Salsburg operator 
—(3IooKoo is a contraction of the corresponding space and therefore —(3~^ ^ o"g(^oo-ft'oo)- 
The same argument applies for the corresponding Kirkwood-Salsburg operator JnK^ as 
well. Q 

Another corollary of Theorem 3-3 can be formulated as follows. Let us define the condi- 
tional, finite- volume free energy density p% as: 

Then 

CoroUciry 3.5 Assume that all hypotheses of Theorem 3-3 hold. Then for any u G Dc{oo) 
the unique thermodynamic limit p'^oiP) — l™ivtooPAr(/3) = liniAr|ooPAr exists, provided 
(3 & H^. Moreover the map H^3 (3 ^ PooiP) is analytic. 

Proof of Theorem 3.3: 

The Banach space is the dual of the Banach space *E^ which is formed from all 
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sequences (';/'„((A, /x, a)„)) of functions defined on {fl^ x Qoo x -R^)®" measurable in the 
charge coordinates a and such they have a finite norm * || ||^ given by: 

* II (^n)n h= fdg\ (A, i^,a)\\ MiK 1^, «)n) | • (3.88) 

n>0 

Define in the space (for suitable ^ to be chosen later on) the following linear operators 

" 1 /■ 

(*i^oo((^m)m))n((A,/X,a)n) =E|ry (3.89) 



/=0 

I 



]l[{(^w[-^^oo{{K,l^i,ai) I (A',/x',a'))] - 1] 



i=l 



tpn+l-lii^^ , At , a ), (Az+1, az+l), ...(An, //n, Oin)) 



and 



n I 

*i^iv((V'm)m))n((A, «)„) = E n / <(^'' «') (3-90) 

' 1 

n[exp[(--£:jv((Aj,//i,Q;j) I (A', //',«'))] - 1] 

'0„_l_i((A', /x', a ), (A/+1, a^+i), ...(An, /^n, ^n)). 

Calculating the corresponding dual operators to the operators *Koo and *Kn we arrive at 
the following equalities 

{*K^f = and (*Xjv)* = X^. (3.91) 

From the estimates 

||i^oolk<r'cxp[CC(/5)] (3.92) 

||i^iv||€<r'exp[eC(/3)] (3.93) 

the continuity of operators *Koo and *K]^ in the space and for any ^ > follows. 

Let us choose a value of ^ such that the family {p%)n C is weakly -* prccompact and let 
further be any of the accumulation points. As a result of the fact that finite-component 
sequences of compactly supported functions in the space *E^ form a dense subspace and 
the local decay properties of £oo (respectively Sn) we easily see that: 

*niv • • *E'^ •*Jn- *njv "s^^-ongly *-^oo • *Joo (3.94) 
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and 



■ *Km ■ *Jn ■ strongly ■ (3.95) 

as A?" t oo and where * H^: *En *-^oo, *In are the corresponding predual of the oper- 
ators Hat, -E'jv Jn whose existence and form of action on the space *E^, can be 
easily determined from the very definitions (3.73), (3. 69), (3. 77) and remark after (3.77). 

Taking into account that 0.% in the weak -* topology of and the convergence 

(3.92) we conclude that any must fulfill the identity: 

P%o = (-/9)/ooi^ooP^ - (3.96) 

Therefore assuming that j3 G H^{(3) we conclude for the existence of a unique thermo- 
dynamic limit which does not depends of cu, provided cu e T>p{oo) and moreover 
Pn ~^ Poo weak -* topology. The same reasoning applies to as well with the 

conclusion that the corresponding thermodynamic limit poo = P^- Standard application 
of the Mayer-MontroU identities (see i.e. [22], [26]) then improves the proven weak -* 
convergence p% — > p^ to the locally uniform componentwise as stated in Theorem. The 
remaining assertions are easier and will not be proven here. 

Effective applications of Theorem 3-3 amounts to the study of the spectral properties of the 
operator Joo-^oo in the space E^. Using the ideas and methods of [27, 28] it can be proved 
that (T^(n(A)^7vi^^n(A) = {-P e C I Z{j{{-f3)-^) = 0}. However the localisation 
and the flow of zeros of the partition function Z{j as A?" f oo is outside the scope of the 
present paper. It is not known to the authors whether the region H{(3) is close to the low 
temperature region. 

It should be mentioned that in the physical literature the limit temperature T f oo and Z | 
oo has been studied before and is known as Vlasov limit in theory of plasma physics [29]. 
The Vlasov limiting region exactly corresponds to the low temperature behaviour of dipole 
systems studied by us. Whether one can extract from the Vlasow theory some relevant 
information to the dipole systems describing fluctuation from the crystalline ground states 
remains to be studied. 
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Chapter 4 



Concluding remarks 



In this section we shall clarify the meaning of the hypotheses that we stated in the course 
of our analysis and put our work in a certain perspective from which further possible 
developments and the meaning of the low temperature expansion introduced here become 
clearer. 

First we shall discuss hypothesis H3. The Taylor remainder (2-41) can be written as foUws 

^a.(A - I ^^) = £ c^tS^pV](A + ty,-i,- ty,) {y, - y,f (4.1) 

Due to the geometry of i?^, even if V has compact support the size of the fluctuations yx 
can be arbitrary. The only restriction is that \ + yx — ii — y^ G supp V . This leads to the 
question whether the Taylor remainder (4.1) can be a bounded function in the variable 
yx — y^- From the elementary estimate 

I A + yA-z^-y^ l^l|A-/^| -\y\-yn\\ 

and the assumption that V is central it follows that this can be realized. As an example 
let us take: 

, , I exp[ 1 12 , o ] for I a; |< £. , ^ 

fe{x) = ^ ,^ ri^l ' ' (4.2) 

10 \ X \> e 

Then from (4.2) we obtain 

1 

I /^(A - n + yx-yu) \< exp[- 



£2- \\X- ii \ - \ yx-y^\ 



if I X- ii + yx-y^i \> e and fe{X- /i + yx-Ufi) = if | \- ix + y\-y^ |> e. This forces an 
exponentially fast decay oi \ fe{X- n + y\- y^) \ \ {yx-y^f \ as \ yx-yfj,\]oo provided 
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A, ijl are kept fixed. 

The partition function (2.30) can be rewritten as: 



(4.3) 

therefore the corresponding free energy density is 



\im ^^/n^T^ = (4.6) 



The results of section 3 enables us to control rigorously 

lim -jr^MZ'^W) ^ f£{P) (4.5) 

N^oo p I Ajv I 

which we call small fluctuations contribution to limjv^ooPjv(/^) = Poo(/^) provided it exists. 
It is reasonable to expect that 

N-^'^ I An I Z%{p) 

therefore the only (but highly nontrivial) problem to solve is to show (with some additional 
hypothesis on [£^^V^]) the existence of the limit 

This hmit be called large fluctuations free energy density. The superstability hke hypoth- 
esis on [D^y] can produce a probability like estimates that make the appearance of very 
large dipoles very unprobable and this might be sufficient for the proof of the existence of 

Let //^ be the corresponding canonical Gibbs ensemble measure for the corresponding real 
system of particles at inverse temperature (5 and with the density p. Assume that the 
lattice C is the ground state configuration. We will say that the local finite configuration 
uj is an almost £-crystalline configuration iff for sufficiently > a map 

(where uj^\ is the restriction of uj to the set A) defined by: 

r^{x,) = r(i) 

where: = A if there exists only one I e C such that \xi — \\— min{| a:;^ — A 1 1 A' e £} 
and if there are several such A e £ on which the minimum of {| Xi — X' \\ A' e £} is 
achieved the A is choosen to be the first of them in the natural lexicographic order given 
to £. 
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exists and is bijective. 



A Borcl subset S of the space of all locally finite configurations is called uniformly 
crystalline iff there exists iV(S) = sup^^^^ < oo. If we suppose that for sufficiently 
small ^ the Gibbs measure fi^ is supported by some uniformly crystalline subset then 
of course the problem of large ffuctuations will be overcomed by a suitable relabeling 
procedura. However at present only a very limited knowledge on the support properties 
of the corresponding Gibbs measures is available, which is not sufficient for solving this 
problem. 

Now we turn to the discussion of HI. In [32,33] Katz and Duneau have showed that 
the set of smooth central potentials with compact support that leads to an energy func- 
tion S{{x)n) possessing Morse function properties is residual in the space of C°°(-R^,-R) 
functions equipped with the Whitney topology. Moreover they showed that the Morse 
character of S{{x)n) is precisely the necessary and sufficient condition for having varia- 
tions of potentials which gives rise to a continuous trajectory of equlibria, preserving the 
crystalline symmetries of the starting equilibrium configurations. The abstract results of 
[32,33] substantiate our motivation for imposing Hypothesis 1 onto our scheme. 

That the set of potentials for which H2 holds is large is expressed in the following propo- 
sition. 

Proposition 4.1 The set of potentials $ e C°°{R^,R) for which H2 is true, is an non- 
void open set in the Whitney topology. 

It would of course be interesting to produce explicit examples satisfying all assumptions 
we made, in particular the uniformity in N involved in HI, H2. 

Some preliminary results on this problem are contained in [21] and a more complete dis- 
cussion of this problem is now in preparation [34] . 

Another important question concerns the stability of the crystalline symmetry in the low 
temperature region in the following sense. In the Pirogov-Sinai theory [35] powerful prob- 
abilistic arguments are worked out to show that the corresponding pure phases do not 
fiuctuate much around the ground state configuration in the low temperature region (This 
is expressed as follows: The probability of the occurrence a large contour F with length | V \ 
is less then exp(— /9 | F |). It is an important problem to develop systematically methods 
which show that the typical configurations of a real system (in the thermodynamic limit) 
in low temperature thermal equilibrium do not fiuctuate much around the corresponding 
crystalline ground state configuration. But an even more important and harder problem is 
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to introduce an equivalent notion of the Peierls stability condition of 'small' perturbations 
of the ground state which leads to stability of crystalline symmetry in the low temperature 
region for the corresponding Gibbs state. 
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